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Using a time-dependent approach the analysis and optimization of a planar FEL-amplifier with 
an axial magnetic field and an irregular waveguide is performed. By applying methods of nonlinear 
dynamics three-dimensional equations of motion and the excitation equation are partly integrated 
in an analytical way. As a result, a self-consistent reduced model of the FEL is built in special 
phase space. The reduced model is the generalization of the Colson-Bonifacio model and takes into 
account the intricate dynamics of electrons in the pump magnetic field and the intramode scattering 
in the irregular waveguide. The reduced model and concepts of evolutionary computation are used 
to find optimal waveguide profiles. The numerical simulation of the original non-simplified model is 
performed to check the effectiveness of found optimal profiles. The FEL parameters are chosen to 
be close to the parameters of the experiment (S. Cheng et al. IEEE Trans. Plasma Sci. 1996, vol. 
24, p. 750), in which a sheet electron beam with the moderate thickness interacts with the TE01 
mode of a rectangular waveguide. The results strongly indicate that one can improve the efficiency 
by a factor of five or six if the FEL operates in the magnetoresonance regime and if the irregular 
waveguide with the optimized profile is used. 
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I. INTRODUCTION 



The recent progress in the theory and experiment of 
free-electron lasers (FELs) and gyrotrons [l], Q with 
Bragg cavities is strongly indicative that the applica- 
tion of novel electrodynamical structures provides the 
opportunity to realize unique properties of FELs to a 
large measure. For example, Bragg cavities, which are 
periodic arrays of varying dielectric or metallic struc- 
tures, stimulate interest in traditional microwave appli- 
cations because they can be built oversized (quasiopti- 
cal) and, therefore, employed at higher frequencies and 
higher power. At the same time the investigation of 
traveling wave tubes (TWT) Q shows that the TWT 
efficiency based on a regular (along the interaction re- 
gion) electrodynamical structure is far from its possible 
maximal value. In fact, the difference between the cold 
phase velocity and the average velocity of the beam is 
the control parameter of the beam-wave interaction. By 
changing this parameter along the interaction region one 
can control the beam bunching and the energy transfer 
between bunches and microwaves. The local variation 
of the cold phase velocity along the region depends on 
the local variation in the waveguide profile. Then, in 
an effort to control the beam-wave interaction and im- 
prove the efficiency one should use irregular electrody- 
namic structures. Specifically, the combination of Bragg 
reflectors Q and the section of an irregular waveguide 
seems to be a promising electrodynamic structure for a 
high-efficiency powerful FEL. The analysis of a planar 
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FEL-amplifier with an axial magnetic field and an irreg- 
ular waveguide is the topic of the present paper. I focus 
my attention on this FEL configuration because it is well 
known that through the use of the axial magnetic field 
one can substantially improve the efficiency, as the cy- 
clotron frequency tends to the undulator frequency Q ■ 

It is worth noting that in vacuum electronic sources 
of coherent radiation the electron beams are far from 
the statistical equilibrium and during their interaction 
with radiation they remain sufficiently nonequilibrium 
because of the large free length [f| 0]. Thus, the effi- 
ciency of the transfer of the electrons' kinetic energy into 
radiation, basically, may be close to 100% (the klystron 
or traveling wave tube are the examples of high-efficiency 
devices) and the challenge is to optimize the beam-wave 
interaction by controlling the most important parame- 
ters. There are several ways to improve the FEL ef- 
ficiency: optimization of electron beam characteristics 
(for example, development of beams with the optimal 
distribution of the axial velocity across a beam when 
the effect of beam finite thickness is relevant), taper- 
ing of the undulator or the axial magnetic field, profiling 
of waveguide/resonantor walls. In particular, the effec- 
tiveness and reliability of the undulator tapering were 
demonstrated theoretically Q and confirmed experi- 
mentally to a great advantage. In the experiment fioj 
the saturated power of 180 MW (6% efficiency) has been 
increased to 1.0 GW (34% efficiency) by optimizing the 
wiggler profile in such a way that the beam-wave reso- 
nance condition remains fulfilled for many electrons, as 
the electrons lose their energy. The numerical simula- 
tion mentioned in (lo| indicates that one can trap about 
75% of the electron beam and reduce its energy by about 
45%. A high effectiveness of the undulator tapering was 
also demonstrated for a FEL with an axial guide mag- 
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nctic field [1JJ (the maximum efficiency was increased 
by almost 700% as compared to the untapered configu- 
ration). At the same time there exist cases where the 
convenient undulator profiling cannot be used or ensure 
desired enhancement. In particular, if the waveguide 
backward mode is used as an electromagnetic undulator, 
then, clearly, one has to change electrodynamic structure 
characteristics to control electromagnetic undulator pa- 
rameters. The optimization of the electromagnetic struc- 
ture also seems to be more efficient than the undulator 
profiling if the effect of the electron beam finite thick- 
ness is relevant. Indeed, given the FEL amplifier oper- 
ates with Group II orbit parameters, a negative masslikc 
effect occurs [12| in which the electrons are axially ac- 
celerated, as they lose energy to the wave. Hence, the 
electrons must be decelerated to maintain the beam- wave 
resonance. This is accomplished by an upward taper in 
the undulator field [TTJ] . At the same time, the electrons 
with different initial transverse positions are exposed to 
the action of different magnitudes of the pump magnetic 
field. As a result, the average velocity of the electron 
depends on its initial transverse position, and different 
beam layers have different dctunings with the wave be- 
cause the average velocity of the electron governs the ini- 
tial detuning. If the undulator field is tapered 'up', then 
the detuning of external layers of a thick beam increases 
and the contribution of these layer to the total efficiency 
shows a certain decrease. In the present paper I demon- 
strate that one can effectively suppress layering and the 
saturation efficiency effect by using the optimized profiled 
waveguide. I believe this technique to be useful for the 
development of powerful thick-beam FELs (for example 
the FEM experiments [l3[ on generation and transport of 
two intense beams have been performed of late: 0.8 MeV 
electron energy, current densities of up to 1.5 kA/cm 2 , 
0.4 x 7.0 cm 2 beam cross sections). 

This paper is structured as follows: in Sec. II the prob- 
lem statement for the planar FEL-amplifier with the ax- 
ial magnetic field and the irregular waveguide is defined. 
The integrals of motion for a test electron in the pump 
magnetic field are constructed in Sec. III. With these in- 
tegrals and the method of nonlinear resonance the FEL 
reduced model is derived in Sec. IV. In the following 
section the principle of the beam-wave control are con- 
sidered and a practical example of the optimized FEL is 
given. The obtained results are discussed in Summary, 
and, finally, the wave interaction in the irregular waveg- 
uide is examined in Appendices A and B. 



II. THE THEORETICAL MODEL 

Let a sheet relativistic electron beam be injected into 
an irregular waveguide located in the external pump mag- 
netic field that consists of the magnetic field of a linearly 
polarized (planar) undulator and a uniform axial mag- 
netic field (see Fig.QJ. The pump magnetic field is given 



by the vector-potential: 

AP( f) = (B u /k u ) cosh(fc u y) cos(k u z) + B\\y. (1) 

Here B\\ is the uniform axial guide field, B u is the mag- 
nitude of the planar undulator field [l4j . k u = 2tt/X u and 
A u are the wave number and the period of the undula- 
tor, respectively. In numerical simulations we model the 
injection of the electron beam into the interaction region 
by allowing the undulator amplitude to increase adiabat- 
ically from zero to a constant level over N u undulator 
periods [l5[ . The unmodulated electron beam enters the 
interaction region, z G [0,L], with average longitudinal 
velocity Vji . The irregular waveguide boundaries are set 
by expressions: x = ±a/2 and y = ±w(z)/2 (a ^> w), 
where w(z) describes the varying distance between two 
wide walls of the waveguide, and w'(0) = w'(L) = 0. Let 
the FEL-amplifier be seeded by the TEoi mode, which 
is resonant (synchronous) with the electron beam, the 
mode frequency and the amplitude at the input into the 
interaction region (z = 0) equaling us and Vq, respec- 
tively. We consider that the interaction region is ideally 
matched to the regular output waveguide at the section 
z = L. 

Since the narrow walls are profiled, y = ±w( z)/2, one 
can apply the local Fourier-series expansion [ 1 €31 ] over y to 
derive the coupled set of equations governing the evolu- 
tion of TE„ m and TM nm modes (subscripts n and m cor- 
respond to the field variation along the wide and narrow 
walls, respectively). Modes with odd TE„ i0 dd, TM„ jOC jd 
and even TE„ jCVCn , TM„. cvcn variations are not coupled 
because of the waveguide symmetry with respect to xz- 
plane. We will hold that A > 2w(z)/3, and the TE 0m 
modes for m > 3 are then evanescent and the scattering 
of the seed TE i mode to those modes will be neglected. 
We also ignore the electron beam mode generation. Un- 
der this assumptions the evolution of the signal TEoi 
mode is governed by the rr-componcnt of vector-potential 



^ ( ^)=Re{n^)^co S (^)e--}. (2) 
Here V(z,t) is the slow-in-time amplitude satisfying the 
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FIG. 1: Sketch of the FEL in the x — cross section. 
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Xdz 2 



2ik z d 

V„ r dt 



}v<„«)~| 



w(z)/2 pt-\-7T 1 04 

dy I dt' j x (r,t') cos( 



w(z)/2 



t — Tr/u 




where 



k 2 z (z) 



w(z) 



w{z) 

w'(z)N 2 r 
2w{z) 



dxx 



. (3) 



1 



71" 

T 



is the wave number squared, v gr {z,ui) = {dk z /dui)~ l is 
the group velocity, c is the speed of light and S — a x b. 
The boundary conditions read 



fdV ., TA i ov\ 



•I T T 1 3F\ 

V OZ V gr Ot / 



2ik z V , 



(4) 



= 0. 



The microscopic current density is given by the follow- 
ing sum over electron trajectories, [la ]: 



the mechanical momentum p by P = p+ (e/c)(A p + A s ). 
The initial conditions for the mechanical momentum and 
coordinates read: 



Px\t=t e = Py\t=t e = 0, p z \t=t e =£V\\/c 2 

x\t=u = %o, y\t=u = Vo, z \t=t c = o, 



(7) 



where £ is the initial energy of the electron entering the 
interaction region at the time t e . The excitation equa- 
tion ^ along with the expression for the current den- 
sity ([5]) and the equations of motion generated by the 
Hamiltonian ([5]) describe the electron-wave interaction 
in the studied FEL in a self-consistent way. In the self- 
consistent models the main mathematical obstructions 
are due to the nonlinear character of the equations of 
motion, and in order to perform the analytical treatment 
of the FEL we first have to integrate the equations of 
motion generated by the Hamiltonian ([5]). For the pur- 
pose of subsequent analysis, it is convenient to rewrite 
the Eq. © as 



H(f,P,t) = JH 2 (f,P) + W(r,P,t), 



(8) 



X b /2 Y b /2 

j(r,t) = -g- J dx Q J dy 



dte p{z;ry,t e ) 
Pz{z;r j_o,t e ) 



where % describes the dynamics of a typical electron in 
the pump magnetic field 



-X b /2 -Y b /2 t-L/V {l 

5[r± - ?x(z;r± ,t e )]S[t - t(z;f± ,t e )], (5) 

where Iq is the beam current at the input into the inter- 
action region; Sb = XbY^ is the cross sectional area of the 
beam; p(z;f±Q,t e ) and r±(z; r*j_o, t e ) are the mechanical 
momentum and the transverse coordinate, respectively; 
t(z;f±Q,t e ) is the arrival time of an electron at the po- 
sition z\ t e and fj_o = ^10(^012/0) are the entrance time 
and the transverse coordinates, which the electron has 
at the input of the interaction region. The sheet elec- 
tron beam is lying from xq = —Xb/2 to xq = Xb/2 and 
from y = —Yb/2 to y = Yb/2 in the x and y directions, 
respectively. Since the relativistic electron-wave interac- 
tion is being studied, the nonradiatcd fields (space-charge 
fields) are supposed to be negligible. The relativistic ef- 
fects result in that the force f m caused by the nonra- 
diated magnetic field partially suppresses the defocusing 
action of the transversal part of the force f^, pot ^ caused 
by the potential part of the nonradiated electric field, the 
axial component of f^ pot ^ being partially compensated by 
the force fi rot ^ caused by the rotational part of the non- 
radiated electric field (see [13] for details). 

The motion of a typical electron within the electron 
beam can be described by the relativistic Hamiltonian 



n = y m 2 e c A + (cP - eAP - eA s ) =m e ~/c 2 . (6) 

Here e and m e are the electron charge and rest mass, 
respectively; the canonical momentum P is related to 



H 



{cP x - e[(B u /k u ) cosh(fc u y) cos(fc„z) + B\\y]} 2 + 
(cP y ) 2 + {cP z f + mlA V2 = m e7o c 2 , (9) 



and the ponderomotive perturbation reads 

W = -2e(cP x -eA p )A s x + (eA s x ) 2 . (10) 

We start our analysis with the integration of the equa- 
tion of motion generated by the unperturbed Hamilto- 
nian ([9]). The above integration is the nontrivial prob- 
lem because the nonlinear dynamical system © is not 
globally integrable [lH and exhibits chaotic behavior if 
the absolute value of the difference between the normal 
undulator and normal cyclotron frequencies is less than 
the betatron frequency [l7| ■ The dynamics of electrons in 
an ideal undulator and an axial magnetic field was stud- 
ied in details in Ref. [l7| using Lindshtedt's perturbation 
method in configuration space. This allowed one to build 
a linear microwave theory and analyze the electron- wave 
interaction in the magnctorcsonance regime. However, 
to build the nonlinear theory we have to study the be- 
haviour of a test electron in the pump magnetic field in 
a more specific way. In the next section we build the ap- 
proximate solution for ([9]) by means of the superconver- 
gence method in action-angle variable space and derives 
the explicit analytical expressions for the region location 
of chaotic dynamics in the parameter space. 
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III. DYNAMICS OF ELECTRONS IN THE 
PUMP MAGNETIC FIELD 

A. Action-angle variables formulation 

In order to apply the perturbation method to the non- 
linear system with the Hamiltonian (|9]) the latter is di- 
vided into a nonperturbed (integrable) part that corre- 
sponds to the electron motion in the axial homogeneous 
magnetic field and a small perturbation caused by the 
undulator magnetic field. Based on the nonperturbed 
system (the undulator field is absent B u = 0) we can 
introduce the action-angle variables: 



Py = - 



sin^ c , 



2I C 

V = C\l -COSPr. 



c 2 P x 



(11) 



uo c £ k u 
the initial conditions take the following form: 



it=t e 

lt=te 



J I „2 q2 c „1,2 



P 2 £cosh z [k u y ]/(2uj c ), C | {= 



t=t c 

t e 



(12) 



where /3|| = Vji/c, w c = \e\cBu/£ and lo u = k u Vu are the 
partial cyclotron and undulator frequencies, respectively; 
e = \/2ujp/uj u is the dimensionlcss perturbation parame- 
ter, u>p = \e\cB\J {-\/2£ ) is the betatron frequency. 

The electron motion described by the Hamiltonian 
Eq. © is characterized by two degrees of freedom, 
namely, by undulator (subscript u) and cyclotron (sub- 
script c) degrees of freedom. The transverse inhomo- 
geneous of the realistic undulator field does not lead 
to the appearance of the additional betatron degree of 
freedom, but only modifies the undulator and cyclotron 
motion. For instance, if ojf},u> c <C <jJ u then the un- 
dulator and cyclotron frequencies of coupled nonlinear 
oscillations are Q u = {w 2 — 3w| cosh 2 ^?^]} 1 ^ 2 and 
fi c = {w 2 +w 2 cosh^/cuj/o]} 1 ^ 2 ■ Note that in purely undu- 
lator field (£?|| = 0) the undulator and betatron degrees 
of freedom become split and the dimension of the dynam- 
ical system is also equal to two [l9[ . Using (|TT|) and some 
algebra we may rewrite ([9]) as 



H= ymlc 4 + (ck u I u ) 2 + 2Lu c £I c + EV{I^). (13) 

Here e V(I, 0) is the nonintegrablc undulator perturba- 
tion 



eV(lJ) = J2Y1 V n , m (I c )(cos[n# u + md c ] + 

n— m— 

cos[n0 u -m0c]), (14) 
where for odd values of m the coefficient V n , m is 
Vo.rn = V 2 , m = s 2 P 2 £ 2 sinh[2A]J m [2B] /4, 
Vi, m = e(i\\£^/2uj c £I c cosh[A] (/ m+1 [B] + I m -i[B]), 



and for even values of m the coefficient V n ^ m is 
V , m = e 2 /3 2 £ 2 (<5 rn , + cosh[2A]I ro [2B](2 - <5 TOj0 ))/8, 



Vo 



v . 



—ef3\\£\/2u} c £I c sinh[A] x 
(I m +i[B] + I m -i[B](l — S m! o)), 

Here, A = -(k u y + (\/2up/u! c ) cosh[k u y )) = const, B = 
ck u y/2I c / (ui c £), S m . n is the Kroncckcr delta, I n {x) is the 
modified Bcsscl function of the first kind of order n. The 
equations of motion read: 



Iu 



EE 

n— m— 



'k^I u /£, 



n \ V n ,m(Ic) 

2£ 



sin[?i0 u + to0 c ]± 
sin[n0„ - m0 c ]), 



(15) 



2£ dl c 



i[nd u - nv& c ]). 



The equation set (fT5j) has a lot of internal nonlinear 
resonances (0 U ps ±m?9 c , 20„ ps ±m0 c ) between the 
undulator and cyclotron degrees of freedom. Actually, 
the successive iterations give in the zero approximation 
/ (0) = 4 =tc and tfL 0) = u u (t - t e ), = 7T + w c (i - t e ). 
One can check that the first approximation leads to 
0*( 1 ),/( 1 ) oc e i( - nu '" ±mu 'J t /(nu) u ±mu) c ). As a result, the 
close is nuj u ±mu) c to zero, the more perturbed dynamics 
is. Applying the nonlinear resonance technique to (|15p 
and analyzing each internal nonlinear resonance sepa- 
rately we can show that I c oc y/e in the vicinity of the 
resonance. Using this estimation we can compare the 
levels of the dominance of different resonances (e <C 1, 
k = 1,2,3...): 

if $ u ps (2k - 1)0 C then e V oc e fe / 2 + 3 / 4 C oshL4], 



if § u ps 2fc0 c 



then e V oc e 



fc/2+1 



sinh[A] 



if 2tf u ps (2fc - 1)0 C then e U oc e fe / 2 + 7 /4 smu [2,4], 

if 2tf u ps 2M C then e U oc e fe/2+2 cosh[2A]. 

(16) 

It will be further shown that if e and yo exceed some 
thresholds, then there exist regions of chaotic dynamics 
of the test electron in the phase space. These regions 
of the phase space correspond to the regions of the non- 
linear resonances between degrees of freedom. The most 
important resonances, in the vicinity of which the on- 
set of the chaos can occur, are $ u 
0„ ps 30 c (see Fig. [HO. 



C , U « 20 c and 



B. Superconvergent method 

An efficient method for analytical treatment of Hamil- 
tonian systems is the application of canonical transfor- 
mations to a Hamiltonian [2(|, [Uj], [22[. Then we seek 
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for canonical transformations to new dynamical variables 
such that a new Hamiltonian %' is a function of ac- 
tion variables only. Therefore, new actions become in- 
tegrals of motion. According to the superconvergence 
method [2^ we choose successive canonical transforma- 
tion {I J) -> (hJi) -> (I2J2) -> ■ ■■(!,$) in such a 



way that every next perturbation becomes the order of 
the square of the preceding one: eV — > e 2 V\ — > e 4 V 2 — >■ 
. . . e 2 V n . After two successive canonical transformations 
the Hamiltonian takes the form: 



u 



1 4 

TO„C 



£ 



2 o2 



: 2 cosh 2 [k u y Q ] 
2(A 2 -d 2 ) 



2A c a 1 



! ( A u 



2A 2 d 2 



d 2 ) cosh 2 [fc u y ] 



4d 2 (A 2 -d 2 ) 



2^,2 



e 2 A 2 sinh 2 [fc„j/ ] 



4c 2 (A 2 



4^) 



Ml 1 / 2 

)]} (17) 



Here A u = k u cI u /(/3u£) and A c = k u cl c / (I3\\£.) are the 
dimensionless integrals of motion with an accuracy 0(e 2 ), 
o"o = u c /uj u . The oscillation frequencies are: 



an 

dl c 



dU 

Bin 



£ 2 sinh 2 [fc M y ]A 2 
4d 2 (A 2 -4d 2 ) 2 



+ 



e 2 cosh 2 [k u y ] 



At + 2A 2 u a 2 - 4 



£ 2 <j 2 cosli^i/p] 
2(A 2 -d 2 ) 2 



4d 2 (A 2 -d 2 ) 2 
2e 2 A c q sinh 2 [fc M yo] 
AS "4a 2 

4AgA e 
d (A 2 -d 2 ) 



(18) 
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The velocity components that are needed in the sequel 
are expressed in terms of the action-angle variables as 



Px 



ep\\A 2 u cosh[fc u y ] 
Al^ 2 

p\\ - ep|| cosh[fc u y 
■ cos[#„ - 



cos 6 U + p|| \/2(jqA c cos 
\/o" A, 



\/ct A c 



V^(A„ - 0-0) 
a i , eA„ cosh[fc M y 



V2(A u + ct ) 



4(A 2 -d 2 ) 



cos[6»„ - 9 C ]- 
cos 20, 



(19) 



Here p\\ = £V\\/c 2 and pii = p\\£l u / (k u V\\) are the initial 
and average axial momenta. We have taken account of 
the first non- vanishing corrections only. To complete the 
study we have to determine the approximate adiabatic 
invariants A u and A c . Using the relations between old 
and new variables via the generating functions and ini- 
tial conditions (fT2"f one obtains the set of equations with 
respect to unknown A u and A c . This set has the bulky 
form, and we did not write it here. Instead of this it is 
convenient to introduce two new auxiliary functions x 
and a such that: 



A„ 

— = 1 

X 



_ 2 a 2 cosh 2 1/0] 
' 2{x 2 ~a 2 ) 2 ' 



A r = 



e 2 >f 4 cosh 2 [fc M yo] 

2cr (x 2 - 



(20) 



where x and a satisfied the set of equations: 

f_ 2/ (2^ - g 2 ) cosh 2 [fc M y ] cosh[2k u y ] 

(Jo V 2(x 2 -<7 2 ) 2 4d 2 

sinh 2 [fc„2/ ]\ 1 £ 2 cosh 2 [fc„?/o](3x 2 + d 2 ) 



x 2 - 4cr 2 



x=l 



4(x 2 ~<t 2 ) 2 



(21) 



Then the frequencies (fT8")) are expressed in terms of un- 
known constants x and a in a simple way: 



(22) 



The average axial momentum and velocity prove to be 
equal p\\ = xp\\ and v\\ — xV\\. Recall that e = \/2ujp/u) u 
and (To = w c /oj u . 

Let us find the approximate solution to equation set 
(|2"TT) . For further analysis let us assume that x ~ 1 and 
a ~ ctq. Then, in case of do <C 1 and do 3> 1 we may 
take x = 1 and do m the right-hand sides of Eqs. (|21|) to 
obtain the explicit solution. To consider the case do ~ 1 
we introduce a new small magnitude \l = x — a, /i <C x. 
Neglecting fi in such expressions (/x + x) we get the cubic 
equation with respect to \x: 



(j? + cin 2 + c 2 /i + c 3 



0, 



(23) 



where c\ = do — 1, c 2 = £ 2 do cosh 2 [fc„yo]/4, C3 = 
£ 2 (2 + d ) cosh 2 [fc„?/o]/8. The discriminant analysis 
£>(£,d ,2/ ) = p s /27 + q 2 /4 (p = -c 2 /3 + c 2 , g = 
2cf/3 3 — C1C2/3 + C3) of cubic equation (|23j) shows 
that D(e,ao,yo) < in the region do < dg™* 1 , and 
D(s, do,jyo) > in the region do > dp™* 1 . The quantity 
do"' 1 (£,yo) is the solution to equation D(s, do,yo) = 
and equals: 



d n =1 



9£ cosh[/c u y ] \ 2/3 



3 /£cosh[fe u y ] \ 4/3 



6 



7£ 2 

— ■ (24) 
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In the region ao < (Jq the solution of equation ([23 
has the following form: 

(j, = -(1 - cr )/3 - 2^J-p/3 cos[(c 



2tt)/3], (25) 



where cos a = — q/[2\J— (p/3) 3 ]. In the region a® > 
ctq 1 " 1 ' 1 the solution of equation ([23]) reads: 



/' 



1 - (T 



3(v^D-g/2) 1 /3 



+ (V^-q/2) 1 / 3 . (26) 



For (Jq < e the above-mentioned e expansion is not quite 
correct and this case should be treated separately. The 
analysis shows that the trajectories remain unchanged 
but to calculate the cyclotron frequency we have to make 



use of another formula O c = {i 



cosh [k u y ]} 



1/2 



The comparison of the results for Cl u and Q c obtained 
by using the analytical expressions (solid lines) and the 
numerical simulation of Eqs. (|15[) (dots) is demonstrated 



1.5 



G 

T3 



3 

3 
G 



§- 



1.0 



0.5 



3 

• 3 


„ 

a 

CO 




.G 


i 

cf 






REGION; 








' 3 


y 




'3^ 


o 

u 





CHAOTIC 
REGION 
AROUND 
NONLINEAR 
RESONANCES 

and 

Q„ ~ 2Q 




S = 0.2 
y = 0.05 



0.25 0.5 0.75 1.0 1.25 1.5 
Normalized Partial Cyclotron Frequency <d c /k> u 

FIG. 2: Undulator and cyclotron frequencies vs. the par- 
tial cyclotron frequency. The adiabatic undulator entrance of 
electrons to the interaction region is neglected. Solid lines are 
for the analytical results, while dots correspond to the results 
of the numerical simulation. 
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FIG. 3: Undulator and cyclotron frequencies vs. the par- 
tial cyclotron frequency. The adiabatic undulator entrance is 
taken into account. 



in Fig. [2] and |3l For both of the situations the analyt- 
ical and numerical results are seemed to be in a good 
agrcment. Note that in our analytical study we ignore 
the adiabatic undulator entrance of electrons to the in- 
teraction region. The comparison between Figs. [2] and [3] 
gives a clear indication that the adiabatic undulator en- 
trance "improve integrability" of the nonlinear dynami- 
cal system (|15[) and reduce the width of chaotic regions. 
In [17| it was found that the initial positive value for the 
x-component of the velocity leads to the suppression of 
the chaotic region. A test electron acquires such a posi- 
tive average correction to v x passing through the region 
of the adiabatic entrance. And, as a result, the chaotic 
dynamics of the electron in the regular undulator region 
is partially suppressed. A rough analytical estimate in- 
dicates that the average correction is about one fourth of 
the amplitude of i^-oscillations in the regular undulator. 



C. Chaotic motion 

The equation set (fT5"|) has a lot of nonlinear reso- 
nances (fTrj|) [n&u ~ tw&c) between the undulator and 
cyclotron degrees of freedom, therefore one can expect 
the appearance of chaotic dynamics in the system be- 
havior. In Fig. 2] we have demonstrated the Poincare 
mapping in which the primary and the higher resonances 
are seen; the separate dots correspond to the stochastic 
trajectories. Note that the average axial velocity for the 
majority of the stochastic trajectories equals zero. As 
evident from Eq. (|2T)|) and Eq. (f2"Tj) . the undulator action 
(integral of motion!) can vanish under some conditions. 
This means the destruction of the integral motion [23|, 
ch. 5] and chaotization of the dynamics of the test elec- 
tron. Hence, the motion becomes stochastic if the dif- 
ference between the undulator and cyclotron frequencies 
becomes less than the betatron frequency 



|fi„ — O c | < V2LO/3 cosh[k u yo] (|« 



<j\ < scosh[k u yo])- 
(27) 
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FIG. 4: Poincare mapping to the system (|15|) ; #c = 0. 
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FIG. 5: The major Lyapunov exponent map. Black solid 
lines show the boundaries of the chaotic region according to 
the analytical formulae (|28l) . 



Such a criteria was initially proposed in [171 ] using some 
numerical findings. With the derived tt and a we get the 
expression describing the location of the region of the 
dynamic chaos: 

critl ^ _ ^- crit2 
<7 < (TO < CT Q , 



^.crit2 



28e co$h[k u yo] be 2 cosh 2 [k u y ] 

(28) 



27 



18 



It turns out that the chaotization condition is inconsis- 
tent with the solution of equation set (|2"Tj) for e less than 
the minimal value of e mm 



coshffc, 



2/o J 



0.0786. 



(29) 



In the used approximations this implies that there is no 
the chaotic region for e < e mm . In Fig. [5] we have il- 
lustrated the results of the numerical calculations for the 
major Lyapunov exponent. The solid lines calculated 
using equation ([28]) show the boundaries of the chaotic 
region. They are in a good agreement with the results of 
numerical simulation. 

In what follows we analyze the studied FEL within the 
region of the regular dynamic. 



IV. THE REDUCED MODEL OF THE FEL 

In this section we construct the time-dependent re- 
duced FEL model that allows for the complicated dy- 
namics of electrons in the pump magnetic field and in- 
tramode scattering in an irregular waveguide. Since we 
are going to analyze the resonant interaction between 



the microwave and electron beam we may approximately 
rewrite vector-potential A s x as (see [B]) 



r / k% 



(Jiy) e -«ot+m*)}. ( 3o) 

\w(z)/ > 



cos | — — ) e 

^w(z) 

Here = /* k z (z')dz', V+(z,t) and V-(z,t) arc the 

amplitudes of the forward and backward waves governed 
by the following excitation equations 



dV+ 
8V- 



1 dV+ 

Vg r (z) dt 

1 dV- 



J e (z,t), 



(31a) 



dz Vg r (z) dt 



2 ^ v - = ttIt- ( 31b ) 

oz 2k z 



The Eq. pia[) describes the interaction between the beam 
and forward wave, whereas the Eq. (|31b[) describes the 
scattering of the forward wave to the backward one. The 
boundary conditions are 



v+\_ n = v , 



I 2=0 

The effective current 



4H uj / k% 1 



V- 



n/2 



o. 



t-t p +7r/cj 



(32) 



k®Sc y k z w Yf, 



dy Q / dt e x 



-Y b /2 t—tP—ir/u 



Px{z\ y , t e ) cQ J Ky(z; y ,t e ) \ aiut (z;yo,t e )-W(z) 



Pz(z;y ,t e ) 



w(z) 



(33) 



at moment t and position z is generated by the group of 
electrons entered into the interaction region during the 
time interval from (t — t p (z; y ) — n/uj) to (t — t p (z; yo) + 
ir/uS). Here t p (z;yo) is the arrival time of the electron, 
which moves in the pump magnetic field ((T|) , to the cross- 
section z. According to Eq. ((3|) and Eq. ([5]) the integrat- 
ing over t' yields a non-zero result only if t — tt/lo < 
t(z;yo,t e ) < t + tt/uj because of the Dirac delta. Since 
the right-hand side of excitation equation © is a slow 
function of time we may write t(z; yo, t e ) « t e + t p (z; yo) 
and find integration limits with respect to t e . 

In the previous Section we have studied the nonlinear 
system with the Hamiltonian @ and found trajectories 

as functions of actions I and angles d. Now let us take 
into account the ponderomotive potential W(r, P, t). We 
will hold H(r,P) as the integrable part of the Hamilto- 
nian ((HJ and consider the relation (fT9|) between (r, P) and 
(I, f9) as a variable replacement rule, regarding (I, t?) as 

new unknown variables. The perturbation W(I, t) pe- 
riodically depends on angles and # c , therefore it can 
be represented as a double Fourier-series over -& u and -d c : 
W oc w^ m e- luJt+i ^^ +m ^^ +imS ^ (m and n arc inte- 
gers). As we have already known from Sec. lIIII thc slower 
is phase changing Lot — ^ — nd u — md u , the stronger is 
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the action of Fourier-component W n>m (see text below 
Eq. (|15[) ). The main principle of the nonlinear resonance 
analysis is simple [23, ch. 3]: the 'troublesome' resonant 
term is separately extracted from the perturbation ex- 
pansion and, in the sequel, the dynamics caused by this 
term is studied. Further on, we analyze the undulator 
resonance 

d u + k z (z')dz' « uit {(k z + k u )v {l ^Lj). (34) 

Jo 

and may write the ponderomotive perturbation as 
W(I, I, t) « Rc{W(I Ul z@ u ),t)eW*<»% 

where 



W 



£W|A 2 cosh[fc M j/ ] 



V+ + V' 



I k% /ny 
- — cos( 

k z w V w 



reduced FEL model in the following manner: 
^ = -Y(y )Z(0 Re{(F + + F_)e^}, 



dF, 



i8 z (C)F + + 

\v„ r 



1 



3F+ 



dC 

- + i(2k z t 9 -6 z )F.-(l + ^' 

\ Vg r 

z(0 



dr 



d( 2k z 
(36) 



dyUd^ Y(y> )e-^->y'o^o\ 



- F , F_ 



<5y(j/o) 



«||(0) v \\(yo)' 



= 0. 



Here £ = and r = (vnt— z)jl g are the dimensionless 
longitudinal coordinate and "retarded time" [24[; 



F+ = - 



ieao 2 £ g ujp 



V±e 



(37) 



Quantity W(I U , z ($u), t) is a slow function of d u and i 
via a slow variation of the waveguide profile and the mi- 
crowave amplitude. The equations of motion is 



1 8W 



£ 



ic = o, s c 



n u {i) 



fic(i) e/u. 



1 dW 

mW u ' (35) 



To obtain the equations of motion in the simplest form 
and clearly demonstrate the physics of FELs with the ax- 
ial magnetic field we additionally suppose that H— £ <C £ 
('Compton limit'). Applying the method of the nonlinear 
resonance [23|, ch. 3] to Eqs. (|35j) . using z as a new inde- 
pendent variable and defining the ponderomotive phase 
as ip = k u z + f* k z (z')dz' — uit, we get the equation for 




d 2 ip 
Ih 2 



dk z 
dz 



1 U> 



1 dW d 
k u JmlhpW u 



H 



where the derivative of the undulator frequency is 



_d_ 



H 



- - a 



U7, 



u u £ d I u 
4° dl u H " f3p 2 < 



= 1 - 



Expanding the ponderomotive current into a series with 
respect to angles and taking into account only the reso- 
nant term in the excitation equation one can write the 



is the normalized field amplitude; "0(0 T ! Z/o? i>o) = 
ip((, r; yo, ipo) — J 8 z (C)dC is the ponderomotive phase; 
tpo and yo are the initial entrance phase and the initial 
transverse displacement of the electron's position from 
the undulator symmetry plane y = 0; 



1 I 2ttw 2 ^} 
^Zk^U 2 



(38) 



where I a = m e c 3 /e f=a —17 kA is the Alfven current 
(recall that e, Iq < 0). The parameter £ g is called the gain 
length [24| (the spatial growth rate of the FEL without 
the axial magnetic field is equal to l~ x for zero detuning). 
The explicit dependence of the reduced FEL model on the 
transverse electron's position and the axial position are 
given by the relations 



v ^u(yo) cosh[fc u y ] 
y = , — ; tttt-, — r- COS 



7T2/0 



w{z) 



k% 



^(yo)-^ 2 (yo) 
The dimensionless detuning parameter is 

Sz (Q = e g (k z (c) + k u -^ 



k z {z)w{z)' 



(39) 



Note that our model includes two detuning parameters: 
<5 2 (C) changing along the interaction region and 5 y (yo) 
changing across the beam. The physical meaning of these 
parameters is discussed in detail below. 

The time-dependent model (f3l)|) allows for the in- 
tricate dynamics of electrons in the pump magnetic 
field ([T]), the effect of the electron beam finite thickness 
and the intramode scattering in the profiled waveguide 
(the intramode scattering acts actually as a feedback). 
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Our model (|36|) is exactly coincident with the Colson- 
Bohifacio model [25|, [26| if a free-space case is employed, 
the axial magnetic field equals zero, the electron beam 
is ideally thin and ultrarelativistic. In the above case, 
the field amplitude F+ depends mainly on detuning pa- 
rameter S z if the initial value of F+ is sufficiently small. 
Then the FEL operates efficiently if \S Z \ < 2 (see Fig. 2 
in [13,[23])- Otherwise, our model ([55)1 is dependent upon 
more than one parameter and incorporates some novel 
effects. It is worth noting that using the method of non- 
linear resonance in the action-angle phase space one can 
reduce any Hamiltonian with the resonant perturbation 
to the so-called Universal Hamiltonian of nonlinear res- 
onance [iH ch. 3]. This means that any resonant beam- 
wave interaction can be described within the framework 
of the reduced model that includes the pendulum-like 
equations of motion of electrons and the excitation equa- 
tion of Colson-Bonifacio type. The equations of motion 
have dimensions of one and a half. In our case, the beam- 
wave energy transfer occurs through the undulator de- 
gree of freedom only, whereas the energy stored in the 
cyclotron degree remains unchanged. 

What is important is that the model (|3l)]) depends 
solely on the average axial velocity via the detuning but 
does not depend on the particular scalar components of 
the initial velocity. This results in that the FEL effi- 
ciency is only dictated by the spread of the average axial 
velocity such that 5v\\ oc Sv z + £5v x +e 2 5v y , where 8vi(t e ) 
(i = x 1 y, z) is the magnitude of the initial velocity spread 
{8vi/v\\ <C I) and Sv^(t e ) is the average axial velocity 
spread. Recall that e = y/2uifj/uj u is a small parameter. 
It is clear that one first needs to minimize the initial axial 
velocity spread. The analysis of the FEL indicates that 
the velocity spread changes the efficiency insignificantly 
if the detuning caused by the spread is much smaller than 
unity [24[ . For the ideally thin beam (Y& — > 0) it yields 
the condition 



(fc° + k u )£ g /ie„ 



y2/3 



< 1, 



(40) 



where jJLjj., is the variance of 8v\\ . It was shown in 28[ (see 
also the results of the numerical simulation in Qjj | that 
essential decreases in the sensitivity of the the efficiency 
to the initial beam spread can be obtained if the undu- 
lator frequency is close to the cyclotron one (multiplier 
y~ 2 / 3 attains its minimal value). 

In general, the ponderomotive potential enhances as 
the undulator frequency tends to the cyclotron one 
(Y(yo) increases). This results in a stronger coupling 
between the wave and electrons. Such an effect referred 
to as the magnetoresonance is well-known in the liter- 
ature Q and the recent detailed study [13] confirmed 
the usefulness of such a regime for a planar FEL con- 
figuration. However, the magnetoresonance effect is not 
so effective when the beam has a finite thickness. Elec- 
trons with the different initial transverse positions, yo, 
undergo the action of the different magnitudes of the 
pump magnetic field ([1}. Then the average velocity of 




-0.04 -0.02 0.02 0.04 
Position of Electron Beam Layers, yJX u 

FIG. 6: The relative efficiency (histogram) and the initial 
detuning (red solid line) of electron beam layers vs. the nor- 
malized transverse displacement of layers from the symmetry 
plane y — 0. 



the electron depends on its initial transverse positions yo 
(this dependence particulary strong near the magnetores- 
onance). At the same time the average velocity of the 
electron of the beam governs the initial 'transverse' de- 
tuning S y (i/o) between the electron and the wave. Hence, 
the value of S y (yo) changes across the beam, and the 
contribution of different electrons to the total efficiency 
might be quite different. To demonstrate this effect we 
simulate Eqs. (f36|) for the parameters close to the exper- 
iment [29j and assume there is additional axial magnetic 
field 20 kG as well. In the simulation we split the beam 
into 21 layers in the transverse cross-section. Each layer 
is also uniformly distributed into 50 macroparticles en- 
tering within one wave period. Recall that the physical 
system under study is homogenous in the x-direction. 
The results are shown in Fig. [51 The internal layer oper- 
ates in the regime of optimal (with respect to efficiency) 



detuning S Z (C) 



(27/4) x / 3 and 5 y = but the 



external layers operate with the non-optimal detuning 
because of the variation in 5 y (illustrated in Fig. [6] as the 
solid line) across the beam. 

It should be anew mentioned that under certain condi- 
tions the integral of motion of a typical electron fails and 
dynamics becomes chaotic. Now we can derive the chao- 
tization condition including the effect of the microwave. 
In the microwave saturation region one can hold the av- 
erage value of the undulator action as an integral of mo- 
tion, then the condition (J„) > has to be fulfiled to 
preserve the regular dynamics of electrons and the valid- 
ity of model (f36|) . In the case of the steady-state regime, 
when the beam is thin and the waveguide is regular, the 
improved chaotization condition is 



7(|fi« \F+\ 



(O u - fi c ) 2 



u> 8k u £„ 



2/3 



< 1. 



(41) 



Here we have considered that the microwave field modi- 
fies the undulator action by the additional quantity 
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Besides, we used the constant of the motion to the 
Colson-Bonifacio model [301 



\F, 



\Fo\ 



[ d(' 4 

It is obvious from Eq. (|41[) that the microwave may cause 
the chaotization of electrons even if the dynamics was 
regular in the pure pump field. 

Another important feature of the model ([55)1 is that 
it takes into account the effect of waveguide profiling. 
This effect exhibits the coupling between the forward and 
backward waves because of intramode scattering as well 
as the dependence of the wave number, k Zl on the axial 
position z through the varying waveguide width, w{z). 
As a result, the detuning 8 Z is also a function of the axial 
position and its control can be used to govern the beam- 
wave interaction. 



V. CONTROL OF THE BEAM- WAVE 
INTERACTION: FEL WITH THE 
OPTIMIZED WAVEGUIDE PROFILE 

Now discuss the physical principle of the control of the 
beam-wave interaction. For simplicity we consider the 
steady-state regime and the thin beam. We also neglect 
the backward wave generation assuming that w(z) is a 
slow function of z. Rewriting the complex amplitudes 
of the wave and ponderomotive current as F + = \F + \e la 
and J e = |J e |e ra * and using Eqs. (|36p we arrive at the 
system: 



d 2 ^ 

d(- 

da 
~d< 



dC 2 - -|^+l cos (a + ^), 



d\F+ 
d( 



6 z (()-\j e \sm(a-u*)/\F + \. 



| J e \ cos(a — u*), 
(42) 



Arg(- / 
Jo 



dipo 



\J e 



-iip 



dipo 



Here ip, \F+ \ and a are the unknown quantities governed 
by the differential equations, and u* and | J e | are given by 
definition. The phase of the current m* defines the posi- 
tion of the bunch center in the system of coordinates mov- 
ing with the velocity of the beam 0, p. 160] (see also HH, 
p. 325]). One can see that if the phase shift between the 
current and the wave, <fi — a — u*, belongs to the interval 
from — 7r/2 to tt/2, then the right-hand side of the equa- 
tion governing wave amplitude l-F+l (second equation in 
the upper line of (|4*2"j0 has a positive sign and the ampli- 
tude itself grows. The phase shift 4> governs the energy 
transfer from the beam to the wave (or v. v.) because the 
local interaction power is dP/d( cx | 1 1 .PLi- 1 cos 0. This 
implies that we can increase the efficiency by controlling 
</> along the interaction region by changing the detuning 
parameter 5 2 (£) in an appropriate way. The idea of such 
an optimization was originally proposed in the TWT the- 
ory [321 ] . Here, for example, we demonstrate the simple 




2 4 6 8 10 
Dimensionless Axial Distance C, 

FIG. 7: The simulation results for the FEL with and without 
the phase shift optimization are demonstrated. 



indirect optimization method [32j. Now assume that the 
phase shift, 0, satisfies the relation: 



a — = 



'opt 



(C) = — 7rsin 2 



2L. 



pr 3 



pr 



(43) 



OL U* — 4>opt(() = 0, 



C > L pr , 



where L pr is the start point of the region with the per- 
manent value of 4>. Then we have to find ip and 
using Eqs. (|42|) and Eq. (|43|) (in the right-hand sides of 
Eqs. (|42p the expression (a — u*) should be replaced by 
4>opt)- Now we can restore the information about the 
waveguide profile using the equation for the detuning pa- 
rameter that follows from (l42ll: 



*.(C) 



du* 

~dC 



d<t>„pt 
dC 



+ (\J e \sm4> opt )/\F 



(44) 



The results from the calculation of the amplitudes of the 
wave l-F+l and the current |J e |, and the phase shift <j> 
are shown in Fig. [7J In this figure we also plotted the 
FEL characteristics for the constant detuning. We can 
see that the wave amplitude can be significantly enhanced 
(efficiency increased several times). However, the demon- 
strated optimization technique is useful only for a slightly 
improved efficiency because the waveguide profiles are to 
be rather complicated from the practical point of view 
in an effort to considerably increase the efficiency. Then 
more elaborated mathematical approaches, which simul- 
taneously allow one to control the practical realizability 
of optimal waveguide profiles, should be used. In this 
paper we apply some type of a genetic algorithm [33| to 
perform the FEL optimization. The principle of evolu- 
tionary optimization is rather simple: we generate a lot 
of waveguide profiles and then perform numerical simula- 
tion of the reduced model (f3"6")l using these profiles. Then 
we choose the best profiles, cross and modify them, and 
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FIG. 8: The FEL efficiency and waveguide width vs. the 
interaction length. The results of the non-simplified model 
simulation are demonstrated. 

perform the simulation again. As a result one can find a 
few best profiles. Finally we must check that these found 
optimal profiles are really useful. We have to simulate 
the non-simplified original model (formulated in Sec. |TTJ) 
using these profiles. Let us consider the results of the 
optimization for a practical example. 

The FEL parameters are chosen to be close to the 
parameters of experiment (29[: 450-kV beam voltage, 
|/o| = 16- A beam current, 1.0 mm x 2.0 cm sheet elec- 
tron beam interacts with the TEoi mode (the field vary- 
ing along the narrow wall) of the 4.5 mm x 4.0 cm rectan- 
gular waveguide. The undulator magnitude increases adi- 
abatically within six periods and the undulator is char- 
acterized by parameters B u = 3.5 kG and \ u = 1.0 cm in 
the regular region. A 1-kW input signal with the 4.0 mm 
wavelength is injected. In our simulation we assume that 
there is the axial magnetic field 20 kG as well. The wave- 
length is slightly different from that in the experiment 
because of the different average axial velocity. In Fig.JSJA. 
the results for the FEL with the axial field but without 
optimization are shown. Using the magnetoresonance ef- 



VI. SUMMARY AND DISCUSSION 



The operation of the planar FEL-amplifier with the ax- 
ial magnetic field and the irregular waveguide is studied. 
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feet we can significantly enhance the efficiency. It was 4% 
efficiency without the axial field in the experiment and 
it is 12% efficiency with the axial magnetic field. How- 
ever, there is a weak interaction between the external 
beam layers and the microwave because different layers 
of the electron beam have different "transverse" detun- 
ing with the wave due to the transverse inhomogeneity of 
the pump magnetic field. Geometric positions of different 
beam layers at the beginning of the interaction region are 
shown inside the dotted ellipse. The black curve is for 
the central layer. Other layers arc displaced with respect 
to the symmetry plane y = 0. 

In Fig. [8j3 the results for the optimized FEL with the 
axial field are presented. Using the waveguide with the 
optimized profile one can double the efficiency so that the 
final efficiency is around 22%. We also see the external 
layers interact with the wave much more effectively in the 
optimized FEL. So, by changing the waveguide profile we 
control beam-wave interaction thus increasing the FEL 
efficiency 

I 

The self-consistent model, which includes the excitation 
equation and the equations of motion along with the ex- 
pressions for the radiated field and the microscopic cur- 
rent density, is formulated. In order to find the parame- 
ters and the waveguide profile that provides the maximal 
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efficiency one has to perform some optimization of the 
FEL. However, the conventional numerical optimization 
methods fail to work because a vast amount of compu- 
tational resources is required. Typically, about several 
thousand equations of motions and the partial differen- 
tial equation for the wave amplitude have to be simu- 
lated. In this paper I propose another approach to the 
problem. The investigation is divided into several stages: 
initially I partly integrate equations of motion and the 
excitation equation in an analytical way using methods 
of nonlinear dynamics. As a result, the universal reduced 
FEL model is derived in special phase space. Then with 
this model and some principles of evolutionary computa- 
tions (genetic algorithms) I perform the numerical opti- 
mization of the waveguide profile. Finally, the simulation 
of the non-simplified original model using the found op- 
timal waveguide profiles is carried out. So, one can come 
closer to understanding of what increase in the efficiency 
can be achieved in practise. 

To derive the reduce model one first can find the inte- 
grals of motion of a test electron in the pump magnetic 
field (JXJ) applying the overconvergent method. The dy- 
namics is completely governed by the two integrals of 
motion corresponding to the undulator and cyclotron 
degrees of freedom. At the same time it is reason- 
able to use other two parameters that completely de- 
fine the dynamics as well: the first governing parameter 
e = y2tup cosh[k u yo]/uj u is the level of nonlinearity and 
the second one <7q = ui c /uj u shows how the system is close 
to the magnetoresonance. The complete description of 
the dynamics is given in terms of these parameters. It is 
well known that the d yna mics in the pump field is chaotic 
for some parameters so the explicit expression ([25]) 
describing the region location of the dynamic chaos in the 
parameter space (e, cto) are derived and the existence do- 
main of integrals of motion is formulated. From the plot 
of the Lyapunov exponent we see that analytical formu- 
lae (|28p give an accurate definitions of the chaotic zone 
boundaries. Note that the afore-mentioned technique can 
also be applied to the pure undulator field (there is no 
the axial field). In this case from Eq. (f2"T)l it follows a 
simple chaotization condition that in terms of paper [l9j 
reads 



a w > f3 zb j b /(V2cosh[k w x b ]). 



(45) 



In Fig. [5] the regular and chaotic regions in the param- 
eter space (k w Xh, a v ,) : according to Eq. (|45j) , are shown. 
Comparing Fig. 2 of [19| and Fig. [HI of the present paper 
we notice that the proposed simple estimation is in a rea- 
sonable agreement with the numerical simulation result 
cited in paper (l9j . 

Using the information about the electron's dynamics 
in the pump magnetic field one can subsequently ana- 
lyze the dynamics of ensemble of electrons in view of 
the ponderomotive wave. In a special coordinate system 
(that moves on the invariant torus surface if there is no 
the ponderomotive wave) one can split degrees of free- 
dom and partially integrate the equations of motion [23j j . 




0.0 0.2 0.4 0.6 0.8 1.0 
Effective Beam Width kjc h 

FIG. 9: The boundaries of regular and chaotic regions ob- 
tained using Eq. (J45J) (cf. with Fig. 2 in pi]). 



As a result, the universal reduced model of the FEL ([36]) 
that incorporates the intricate dynamics of electrons and 
the intramode scattering is derived. 

What is important is that there are two types of detun- 
ing in reduce model ([36| : 'axial' detuning that changes 
along the interaction region via the profiled waveguide 
width and the 'transverse' detuning that changes across 
the beam because the pump magnetic field is inhomoge- 
neous and the average velocity of the electron depends on 
its initial transverse position. The "transverse" detuning 
causes the thick beam layering and the degradation of 
the external layers' contribution into the total efficiency 
(the degradation is particulary strong near the magne- 
toresonance). In the present paper I demonstrate that 
this problem and the saturation effect can be overcome 
by the control of the beam-wave interaction. The physi- 
cal mechanism of such a control is that by changing the 
waveguide profile one control the axial detuning and thus 
regulate the phase shift between the ponderomotive wave 
and current. This phase shift defines the transfer of the 
energy between the beam and the wave and its regulation 
allows one to optimize the interaction. 

The practical example of optimization of the FEL, 
whose parameters are close to those of the experi- 
ment [23| , is demonstrated. The simulation results based 
on the non-simplified model (see Sec. II) strongly in- 
dicate that combining the magnetoresonance effect with 
the optimized profile waveguide one can enhance the FEL 
efficiency by a factor of five or six. The efficiency in the 
experiment [29| was around 4%. Applying the axial mag- 
netic field the efficiency has been increased up to nearly 
12%, but about 30% of electrons do not interact with the 
wave because of the initial transverse detuning. Follow- 
ing the waveguide optimization the efficiency has reached 
22%, in particular, due to a much more effective interac- 
tion between the external beam layers and the wave. 
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Appendix A: Time-dependent excitation equation of 
an irregular waveguide 

The evolution of the resonant (synchronous) TEoi 
mode is governed by the x-component of the vector- 
potential A x , which satisfies the wave equation 



1 & 



-in 



We seek a solution to the equation of the form: 



(Al) 



A x (f,t) =Re / V { z,u')J—^cos(^) e^'W. 
Jo V w ^ z > Kw ( z ) J 

(A2) 

Substituting flM]) into ||XlJ) we derive the excitation 
equation for the Fourier amplitude V(Z,uj') 



{J^W)}%o/) = -i v w(z)J _ a/2 



8?r 



a/2 



w/2 

dyj x (r,u') cos 

-w/2 \W{Z) 



dxx 



(A3) 



where k 2 z (z) = (uj/c) 2 - (tt/w) 2 - (w'/2w) 2 (l + tt 2 /3) 
and j x {r,u') = tt" 1 j x (r, t)e %u>t dt. Here j x (r,oj') is 
the Fourier amplitude of the current density. We will 
consider that at the section z = the FEL-amplifier is 
seeded by the TE 01 mode with a frequency of uj and am- 
plitude Vo , and the interaction region is ideally matched 
to a regular output waveguide at the section z — L 



2=0 



z=L 



2ik z V S[uj' - uj], 



0. 



(A4) 



The conditions for the waveguide profile at the ends of 
the interaction region have the following form: w'(0) = 
w'(L) = 0. We assume that j x (f, uj') is the narrow-band 
signal with a fundamental frequency of uj. This means 
that the current density can be written as j x (f,t) — 
Re{ J x (f, t)e~ tut }, where J x (r,t) is a slow function of 
time such that 

Jx(f,t) « - / j x {r,t)e lujt dt* 

j x {r,uj + Auj)e~ ,Au)t d{Auj), Auj=uj'-uj. 
Expanding k z (z,oj') into Taylor's series over uj up to the 



linear term, multiplying the Eq. (|A3[) by e~ lAuJt and in- 
tegrating it over Aw from — oo to oo we derive the time- 
dependent excitation equation ([3]) for the slow in time 



amplitude V(z, t) = j°° oo V{z,uj+Auj)e- %Aut d(Auj). The 
solution (|A2j) and the boundary conditions (|A4[) can be 
rewritten as Eq. @ and Eq. (UJ), respectively. 
Appendix B: Intramode scattering in an irregular 
waveguide 



The Eq. ([3]) is used for the numerical simulation of non- 
averaged model of the FEL, but for the analytical study 
we have to rederive the excitation equation in a differ- 
ent form. We seek a solution to Eq. (|A3|) for microwave 
Fourier amplitude V of the form: 



/ \ / h° \ 1/2 f 



k z (z )dz , 
(Bl) 

(k z = k z (Q)) and we allow amplitudes V± to be functions 
of the axial position z. Applying the method of varia- 
tion of constants we derive first-order equations for new 
unknown functions V±(z) 



dV+ 
dz 

dV_ 
dz 



iRHSe-™ dk z V-e 



2i* 



2^kjk z dz 2k z 
iRHSe™ dk z V + e 2 ^ 



dz 2k z 



Here RHS(z) is the right-hand side of Eq. 
boundary conditions become 



V 6[uj' 



UJ 



V- 



(B2a) 
(B2b) 

U) and the 
(B3) 



Note that Eqs. (TBI]). (|B2|) . (|B3|) formally define the ex- 
act solution to Eq. (|A3|) . Now restrict ourself to the 
case of the resonant interaction of a beam with a for- 
ward microwave. The first term in the right-hand side 
of (|B2a[) describes the above-mentioned resonant inter- 
action and should be taken into account, whereas the 
first term in the right-hand side of (|B2b[) is nonreso- 
nant and might be omitted. Wc will also assume that 
k~ 1 d z k z < 1 and neglect the second term in the right- 
hand side of (|B2ap because it describes the second- 
order scattering effect (according to the boundary con- 
ditions for a source- free regular waveguide V-(z) = 0), 
but we will keep the second term in (|B2bl) . Using 
some algebra and performing the inverse Fourier trans- 
formation we derive Eqs. (|30[) and (|3"Tj) for V+(z,t) — 
V+(z,uj + Auj)e-^ t+l ^ 9 ^ d{Auj) and V-{z,t) = 
£° V^(z,uj + Aw)e- iAwi - iAu8 «* 'd(Auj). 
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